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A special version of multi-dimensional simple waves given in [G. Boillat, J. Math. Phys. 
11, 1482-3 (1970)] and [G.M. Webb, R. Ratkiewicz, M. Brio and G.P. Zank, J. Plasma 
Phys. 59, 417-460 (1998)] is employed for fully relativistic fluid and plasma flows. Three 
essential modes: vortex, entropy and sound modes are derived where each of them is different 
from its nonrelativistic analogue. Vortex and entropy modes are formally solved in both the 
laboratory frame and the wave frame (co-moving with the wave front) while the sound mode 
is formally solved only in the wave frame at ultra-relativistic temperatures. In addition, the 
surface which is the boundary between the permitted and forbidden regions of the solution 
is introduced and determined. Finally a symmetry analysis is performed for the vortex 
mode equation up to both point and contact transformations. Fundamental invariants and 
a form of general solutions of point transformations along with some specific examples are 
also derived. 
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1 Introduction 

Investigation of nonlinear phenomena appearing in a very wide area of pure and applied sciences 
has met extremely extensive progresses and developments. These studies which split into nu- 
merical and analytical considerations are essentially and in most cases related to some nonlinear 
differential equations. In spite of such a huge amount of improvements almost all of these equa- 
tions are still far from being understood well. Among those nonlinear systems, a few interesting 
open problems concern the hydrodynamic type of equations governing fluid motions. Especially 
the Euler and Navier-Stokes equations which reveal a mysterious behavior are being intensively 
studied in two main considerations: The incompressible motion mostly dealing with vortex dy- 
namics and the compressible flow concerning the appearance of discontinuities shocks. Both of 
these problems are somehow related to the debate of the regularity of solutions. In the former 
consideration the aim is to understand the mechanism of occurrence of vortex singularities while 
in the latter the shock convergence is the most important question [TJ [2]. The core of these sub- 
jects is the following principal open question: Starting from an initially smooth flow how can we 
predict the appearance of any kind of discontinuity or singularity in later times? In other words, 
what is "hidden" in the smooth initial conditions which causes the non-smoothness in future? 
Numerical evidences as well as some analytical investigations in special cases highly confirm 
the existence of such hidden facts. It seems that these facts are partly related to topological 
properties and partly to the measure theoretical aspects. 
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The wide variety of the application of fluid motions causes to deal with all kinds of differential 
equations namely, elliptic, parabolic and hyperbolic equations. While many systematic treat- 
ments have been found and developed for parabolic and elliptic equations, hyperbolic equations 
are still out of the frame of any well-defined method. These equations posses some character- 
istic curves or surfaces which naturally have the capability of forming discontinuities. This of 
course makes the nature of the solution to be "local" which means that the continuous solution 
may exist only in some part of the space and in some intervals of time. Generally there are 
many topological, geometrical and analytical unknown features determining the validity of any 
solution which are very difficult to discover under the present human knowledge and thus new 
tools are needed. 

We believe the best way to obtain some result is to study useful special cases guiding us to 
more general statements. An excellent and rich class of solutions for compressible ideal flows 
governed by hyperbolic equations lies in the framework of simple (Riemann) waves which have 
the capability of shock formation [JJ [2] and blow up occurrence [3J. Simple waves constructed 
on the basis of characteristics are clearly local which were discovered first by Riemann in the 
1-D form [3] -[TO] and are still the best analytical tools to achieve shock waves[U [21 El El (III 
I12j . By imposing more restrictions and limitations on these solutions it was possible to build 
multidimensional simple waves [13]- [18 j. Even a more generalization yields double waves and 
multi-waves which give more advanced solutions with more intensive blow up[18j-[25j. 

It is obvious that taking into account relativistic effects highly increases the coupling and 
so the nonlinearity of fluid and plasma motions. Relativistic flows have been known for a long 
time|27]-|29] and especially they are important is astrophysical and cosmological phenomena. In 
addition, under recent technical progresses in laser-plasma interaction, plasma accelerators and 
fusion plasmas, the access to relativistic effects in the laboratory is now very easy. Therefore 
a great attention has been paid to analyze relativistic flows. Again the study of simple waves 
plays a very fundamental role as almost the only available nonstationary exact solution with the 
ability of discontinuity formation. 

A very excellent and complete mathematical discussion on one dimensional relativistic MHD 
simple waves has been reviewed by Shikin[30]. Some solutions of these 1-D simple waves are 
found in many papers. Although a relativistic 2-D double wave solution solution has been given 
only for ultra-relativistic fluids [25] but still we observe the missing of a multidimensional simple 
wave for a fully relativistic flow. This task is the aim of the present note in which the approach 
of Ref. [17] is employed and generalized. 

Physically it is a valid question that why we should consider relativistic fluids while usually 
at so high temperatures matter is found in the plasma form and so one has to take into account 
electromagnetic fields leading to MHD equations. However this is not always true because 
sometimes we deal with neutral fluids like neutron stars. Moreover in the absence of any external 
magnetic field and when the typical length and time for the non-neutrality of the plasma are 
sufficiently less than the length and time for macroscopic motions, the plasma can be considered 
as a neutral fluid with at a very high accuracy. Hence it has sense to consider the ideal relativistic 
flow here. 

In the next section after a brief derivation of relativistic ideal fluid equations, a multidi- 
mensional simple wave ansatz is substituted into these equations and various modes and phase 
velocities relative to the laboratory (fixed) frame are found. In Sec. 3 some solutions for the 
vortex and entropy modes are given only in the laboratory frame. The presented solutions are 
very general and formal and a detailed solution is very difficult and needs to determine the 
initial and boundary conditions precisely. Thus, our solutions are very general including many 
arbitrary functions. In Sec. 4 the equations are rewritten in the wave frame and again some 
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simple typical solutions are given for the all three modes vortex, entropy and sound. Especially 
for the sound since its equations are very complicated in the laboratory frame, it is seen in 
Sec. 4 that these equations in the wave frame at ultra-relativistic temperatures are simplified 
and it will be possible to obtain some formal solutions for it. In sec. 5 we investigate symmetry 
properties and their related topics for the vortex mode equation as a sample equation appearing 
in our problem. Finally a summary and concluding remarks are given in Sec. 6. 



2 Multidimensional simple wave formulation 

Relativistic effects in continuum matters in two aspects: Large macroscopic (fluid) velocities 
and relativistic temperatures at which the mean thermal energy of particles are comparable 
with their rest energy. Both of these aspects are included in the energy-momentum tensor 

Tf = wu k Ui-Pb\ , (i, k = 0,1,2,3) (1) 

where u 3 = (7, tv/c) is the contra-variant 4-velocity and thus Uj = (7, — 7v/c) is the co-variant 
4- velocity and w = e + P in which P is the fluid pressure and e is the internal energy (including 
the rest energy) per unit proper volume (unit volume in the inertial frame in which the fluid is 
momentarily at rest). Therefore w is the enthalpy per unit proper volume. Also v is the fluid 
velocity and 7 = (1 — -u 2 /c 2 )~ 1//2 where c is the speed of light. 
Basic equations consist of continuity equation 

(9 1 d 

— (W)=0, or - — ( 7 n) + V • (n 7 v) = , (2) 

and the vanishing 4-divergence of the energy-momentum tensor 

9 : T l k = 0, (z = 0,1,2,3) (3) 



dx k 1 

where n is the number density of fluid particles in the proper frame. By virtue of thermodynamic 
identity TdS = d(w/n) — dP/n (T is the fluid temperature and S is the entropy per unit particle) 
one can combine Eqs. ([2]) and ([3]) to obtain [30] 



— = _ + v-V5 = 0. 4) 

dt dt w 

This equation can be alternatively considered in place of the zeroth component (i = 0) of Eq. ([3]). 
Thus our set of equations consists of five equations (J2]), (jl]) and the spatial components of ([3]). 
This system of course needs a thermodynamical state equation P = P(S,w). However it is 
found that this system of equations takes a more appropriate form by the use of the following 
useful transformation [30J 

W -u* , p=^-, (5) 



mnc w 

where m is the "mean" rest mass of all particles in the fluid. This transformation makes our 
final system of equations to 

~(p«o) + V • 05k) = , (6) 

_£_ + (K .V). = -,VP. (7) 



3 



^ + k-VS = 0, (8) 

c ot 

P = P(S,p). (9) 
Here k 1 = (re , n) and thus rej = (reo, — «0 and 



k° = ko = yK? + w/p , (re = |k|) (10) 

which follows from the identity u l Ui = 1. 

The special feature of a simple wave in any quasi-linear hyperbolic system of equations is 
that all quantities are considered as functions of only one variable which we call it the phase 
and denote by ip = <p(r, t). In our problem we write 

U = I%), (11) 

where 

U = {ft, K!,K2,K 3 ,S) (12) 

is the state vector of the system. Boillat[15] [T7j showed that for a simple wave it is necessary to 
have 

n = n(cp) , - ,_ , = A = \{tp) . (13) 



\V<p\ ^' ' |Vp| 

In other words, the unit vector n normal to the wave front must be only a function of ip and the 
same is true for the phase velocity A. Finally condition (I13p implies that ip must satisfy [15| [T7] 

G(<p,r,t) = f(<p)+\(<p)t-r-n(<p), (14) 

in which / is an arbitrary differentiable function to be fixed through initial conditions. This 
equation clearly means that level surfaces of ip are flat planes. The functional form of n(tp) can 
not be determined from any equation and so it remains arbitrary to be flexible to fit with a 
given condition. 

There are two noticeable points about these multidimensional simple waves. The first is the 
wave breaking at which the time and spatial derivatives of (p and so all variables diverge when 
F — > provided that 

d<P -%?) n(ip) 

„ dG dfUp) d\(ip) da(v) 1 , N 

dip dip dip dip |Vy| 

Equations (fTo]) are easily derived by implicit time and space differentiations of (fT4"|) . Thus our 
simple wave solution is valid only when F > and at any time and point where F = the 
solution is not correct. The second point arises from the dependence of n on ip which implies 
that for two different values cpi and cp% of ip generally n(y?i) and 11(922) are not parallel and thus 
they have an intersection on a line at which the solution is multi-valued which is not accepted. 
Hence, the domain of the valid solution must not contain such intersections. Both of these points 
demonstrate the "local" character of simple waves. 

For a unidirectional 1-D simple wave where n is a constant vector it is possible for each value 
of ip to calculate the time of wave breaking (F = 0) as t c (ip) = — (df / dip) / dX / dip and the earliest 
time of the wave breaking is obtained by solving the equation (dt c /dcp) = 0[7j. Unfortunately 
such a nice situation does not hold in the multidimensional case when n = n(cp). Let us see this 
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in a quantitative way. Singular points (wave breaking) not only must satisfy the simple wave 
condition flH| but also they should fulfil 



F = 



(17) 



Thus, the wave breaking occurs on the line of intersection of the two perpendicular planes 
G = and F = 0. This line is exactly the rotation axis of the wave front at p when p has an 
infinitesimal growth to <p + Sip. This will be easily seen if we observe that the wave front for 
ip + dip must satisfy 



G(ip + dp,r,t) = 



or G(tp, r, t) + F5p = 



which again yields Eqs. (j!4[) and (I17p . We may therefore conclude (without a rigorous proof) 
that if the wave breaking (singularity) line lies out of the region of the solution, the line of 
multi-valuedness will also lay in that region. Besides, since a line of singularity for each value of 
(p exists at each instant of time, it has no sense to speak about t c (<p). However if the fluid fills 
the whole space HL 3 we can obtain a moving surface constructed at any time exactly from all of 
these singular lines at that time. This surface is in fact the boundary between the forbidden 
and permitted regions relative to a simple wave solution. 

Now we substitute the simple wave ansatz (fTT]) into Eqs. ©-([81) supplemented by Eq. Q 
and then divide each equation by | Vy> | and use (fT3l) to obtain the system of five quasi-linear 
coupled equations 

dip 

where A is the 5x5 matrix with the following elements 



(18) 



A 



a n2 



a no, 
P 




7 K 









:K0 






A , 







_ A n dK 

cP dS 
Psn-l 
P 

-PS"2 
P 

Psn-2 
P x 



(19) 



In the above matrix we have used the following notations 



K, ■ 



i=l 



dP 
dp / s 



(20) 



Moreover, in the calculation of and we must assume w = w(p, S) and use Eq. (Tl0|) to 
express kq explicitly as a function of all five variables U = (p, n, S). 
Equation (|18p has a nontrivial solution only when 

det(A) = , 

which constructs a fifth order equation for A with a triple root 

Ai = A2 = A3 



(21) 



c — = v n , 

K 



(22) 



while the fourth and fifth roots A4 and A5 are the larger and smaller roots of the following 
quadratic equation respectively. 



The triplet root is the phase velocity for the two vortex modes and one entropy to be discussed in 
the next section. The roots A4 and A5 are the phase velocities for the forward and backward sound 
modes respectively. Although these modes have nonrelativistic analogue but they significantly 
differ from the nonrelativistic case. 

Substitution of each value of the phase velocity into (|18p yields some ordinary differential 
equations for U(y>) to be solved. For the entropy and vortex modes these equations are not 
difficult and some formal solutions both in the laboratory frame and the wave frame will be 
presented in Sections 3 and 4 respectively. Since the equations for the sound waves are compli- 
cated in the laboratory frame, we go to the wave frame but still they are difficult to solve and 
finally when we consider the physically common case of ultra-relativistic temperatures it will be 
possible to obtain some solutions in Sec. 4. 



3 Vortex and entropy modes 

If we substitute the triplet root A = c^- into (Tl8|) we obtain 



+ ( n " — K ) • — " — UF— = » ( 24 ) 



kq dp dip \ Kq J dip kq dS dtp 



2 dp ( n dS dP 



a ~r + p s~r = 1- = => p (f) = const > ( 25 ) 

dip dip d<p 

0~=0. (26) 
dip 



In Eq. (I25p we have used the second and third equations of (I20|) together with ([9]). Equation 
(|26p admits the two cases of constant entropy (vortex modes) and variable entropy (entropy 
mode). 



3.1 Vortex modes 

We have dS = or 



S(<p) = const , (27) 



which together with (|25p and the state equation ([9]) yields the constancy of p and so all thermo- 
dynamical variables. Thus, only the fluid velocity considered in k and n change with <p which 
n(</?) is an arbitrary suitable function. Regarding the above results in Eqs. (|24l) and ()10p one 
can obtain the equation for K,(ip). 



Q _Kn K \dK =Q ^ (2g) 



Kq J dip 



kq = a/k 2 + wq/p q , (29) 



which must be supplemented by 



where wq and po are constant throughout the wave. The factor — ^t^J in (|28p can not be 

zero because if it is zero we can make its inner product with n and obtain k\ = k 2 = Kq which 
is impossible by (|29p . Therefore Eq. (|28p is equivalent to 

^ = XMx(„-| K ), (30) 
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where X(y>) is an arbitrary continuous function. It is possible to choose two functions Xi((/?) 
and X2(</?) where Xi(^) • X2(<£>) = which gives two perpendicular and independent vortex 
modes similar to the nonrelativistic case [17]. 

It is also worth noting that we can define a generalized vortex 

ft = V x k. = Vw x — = Vdn x — , (31) 

dip dip 

which is constant not only on the wave front but also in advection with the fluid velocity: 

_ + ( y .V )S! = |V v |(^-Aj^ = 0, (32) 
This fact is consistent with the "frozen in" condition of filed lines of fi[31] 

in which the number density in the laboratory frame 777, is constant because 



V • v = cV • — = n - — S-k • — = , 34 



according to Eq. (|28jl . Finally by (]3ip we find 

n v = iv^in • n— = o , 

dip 

by which Eq. (f3"3"J) reduces to 

Equation (130)) has many solutions since X(</?) is an arbitrary continuous function. It is 
therefore easy to choose some suitable simple forms for X such that Eq. (|3U|) can be easily 
solved. As an example consider two different forms 

K 2 

X = an , or alternatively X = a— n , (35) 



K 



n 



where a is a dimensionless constant. Selecting each value form for X from (I35p causes to reduce 
Eq. flUD to 

— = aK x n(ip) . (36) 

dip 

Regardless of the functional form of n(^) it is obvious from (|36|) that | | is constant and 
from (|29p kq is also constant and only the direction of k(c/?) changes by (p. To obtain a more 
special solution let us consider a 2D simple wave with[17j 

n(ip) = (— sin ip, cos cp, 0) , (37) 

in a Cartesian coordinate system. In this case we have 

K = K„n + K t t + K 3 Z , (38) 

where 

K n = k ■ n = — «i sin ip + K2 cos ip , Kt = k ■ t = — (ki cos p + K2 sin ip) , (39) 
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in which t = (— cos p, — sin p, 0) is normal to n. Substitution of (J38J) into (|36J) and using the 
relations 

dn dt 

— = t , — = — n , n x t = z , (40) 

dip dtp 



one finds 



^ = K " < 41 > 

dut , . \ 

-± + K n -aK 3 = , 42 
dtp 

dK3 (A1\ 

— — = -ctK t . (43) 
dp 



Equations flUJ) and flU yield 



k 3 = -«K n + k , (44) 
where R is a constant with the dimension of k (velocity). Then we substitute Kt from (|4ip and 



+ (1 + cr)^ - a/c = , (45) 



k 3 from (jMj) into Q42J) to obtain 



d^ 2 

with a general solution 

K n = K n cos [VI + « 2 (p + 0)] + t~; — 9^' ( 46 ) 

1 + cr 

where K n is a constant with the dimension of k (velocity) while (5 is a dimensionless constant. 
Then from (|4ip we have 



«t = -\/l + a 2 R n svo[y/l + a 2 {(p + /?)] , (47) 

and from (JMJ) we find 

k 3 = -aR n cos[\/l + a 2 (y? + /?)] + 1 k, (48) 

1 + Q z 

It is then straightforward to find K\ and K2 by the use of 

K% = — (n n sin p + ttt cos , K2 = K n cos — Kt sin 93 . (49) 

As mentioned before a complete solution needs more detailed informations about initial and 
boundary conditions which is not of our interest here. 

3.2 Entropy modes 



Here we have dS 7^ in (J26T) and thus p is not constant although by (j25j) P is still constant. It 
is possible to reduce Eq. (J24"P to 



which gives 



n- =0, (50) 

1 ayj ay? ' 



dK _ r , . d(ln«o) 
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where Y(y>) is again an arbitrary continuous function. Let us again choose a suitable form for 
Y(c^) to simplify the solution. For example if Y is parallel to n we see from (1511) that 

k = k C , (52) 

where C = (01,02,03) is a dimensionless constant vector. Since P is constant, the enthalpy w 
becomes only a function of p and thus Eqs. (j52j) and (fTUj) yield 



, c , (53) 

which is valid if |C| < 1. It remains to specify the dependence of p on cp. This dependence 
is arbitrary because the fixing of <p is under our control and we can assume it as an arbitrary 
function of one or more physical variables [17j. 

4 Simple waves presented in the wave frame 

Sometimes apparent forms of those equations concerning a simple wave solution are so compli- 
cated and not soluble easily. A mathematical trick here is to rewrite all equations in terms of 
physical variables as measured in the wave frame. A wave frame depends on the special value 
of the phase ip that is, for any value of p there is a plane wave front defined by Eq. (|14|) moving 
with the phase velocity \ p h = X(p)n(p) and we consider a Lorentz transformation from the 
laboratory frame to the frame co-moving with this wave front. It is thus clear that to each value 
of (p corresponds a unique wave frame. 

We denote all quantities in the wave frame by a prime except scalar quantities such as p,w,n 
etc which are either Lorentz invariant or defined in the proper frame co-moving with the fluid. 
Therefore for the 4-vector k 1 we have 

kq = cosh^KQ + sinhf;fi4 , 

K n = cosh£fi4 + sinh^KQ , (54) 

K± = , 

in which n' n = n' ■ n , = k! — n' n ii and £ depends on p thorough 

tanh e = ^M. (55) 

c 

In the following subsections we apply this method for the vortex, entropy and sound modes and 
give simple formal solutions for each one. 

4.1 Vortex modes in the wave frame 

For this mode we already had 7 = ^ which by the substitution from (I54|) and (|55p into it we 
find 

< = . (56) 

It is therefore seen that in the wave frame the phase velocity takes the simple form through 
the above equation. Then we substitute (f54"|) and (f55j) into ([28]) and use ([56]) and the identity 
k,' ± • n = to obtain 

, 9 . dn\ sinhf . dn', . ^dn'r, . , _ df „ . 

cosh 2 £n • — ± -^k'i • — ± + smh£— ■ ^ + k' cosh£-^ = . (57) 

dp> Kq dp dp dp 
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Since (kq, k) is a 4- vector, Eq. (|29|) is invariant and due to (|56|)we have 

k' = yj k'I + w /po , (58) 

by which Eq. ([57]) reduces to 

«r? = -^ = 47r- (59) 
dip d(p cosh £ ai/? 

Assuming n((/j) is a known function, k' ± and £ should satisfy Eq. (|59|) which obviously admits a 
large amount of freedom. As a very simple solution let us assume the restriction 

%—4n, (60) 

which says that k'^ and k'q are both constant and thus Eq. (|59p can be easily solved to give 
rt 



/ — — = arctan(sinh£) — arctan(sinh£o) = -^riv ~ Po) • (61) 

Jen coshf k 



From the above solution one can find £ in terms of p by which through Eq. (|55f) we have A(</j). 
For k'j_ we formally solve Eq. (fBTIj) : 



' ± f n(<p'W + *'±(<Po) ■ (62) 



Here we gave a very restricted solution just as an example to show the procedure of the solution. 
Depending on initial and boundary conditions it is in principle possible to find more realistic 
solutions although it seems to be very difficult. 



4.2 Entropy mode in the wave frame 

Since the phase velocity £ = ^ is the same as for the vortex mode, Eq. (|56jl is again valid here 
and transforming Eq. (|50p in a manner similar to that performed for the vortex mode we again 
obtain Eq. (|59|) for the entropy wave too but here since P is constant, w is only a function of p 
and thus 

k' = ^K'l + w(p)/p . (63) 

Hence, the meaning of (I59p for the entropy mode is different from this equation for the vortex 
mode. Again as a very restricted simple solution we suggest Eqs. (f6"U|) and (f6"2"|) for K' ± (ip) 
resulting in the constancy of k' ± and assume a given form for p(<p) by which from (j63H we have 
k' (<p) as a known function of ip and thus Eq. (|59p has the formal solution 

arctan(sinh^) — arctan(sinh^o) = ^'i / f ■ (64) 

J<po W') 

4.3 Sound mode in the wave frame 

At first we see that in the laboratory reference frame there are five equations included in (|18|) 
but due to Eq. ([2T]) we have only four independent equations. The last (fifth) equation of (fT8|) 
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by the use of (|19p implies Eq. (|27|) also valid for the sound mode which its substitution into the 
first four equations of (fTHj) yields 

A \ dp \„dK dn 

i k - p— + pn ■ — = , (65) 

c J dip c dip dip 



as the continuity equation and 



A \ dn a 2 dp , . 

-Ko)-i- + —-fn = 0, 66 

c J dp p dip 



as the momentum equation in which a 2 is defined from (|20|) . Equations (|65p and (|66p are four 
equations but only three of them are independent while the phase velocity is A = A4 or A = A5 
which are the roots of the quadratic equation (j23j) . 

Let us rewrite Eqs. (|65|) and (|66|) in terms of the wave frame quantities through Eqs. (|54|) to 
find 

, dp ~ 1 dtl ~dn' n _ cZft'x 



^ + ^ + ^ + coshe ^ n • ~dV = ' (67) 



and 



, . a dp dn j_ d . , . . . dn 

cosnf — r-^n H h n— - (k„ cosh 5 + K n smnf) + (n„ cosnt + K n smnt )— — = , 

pn' n dip dp dp u y v n u (iy? 

with only three independent equations. 

It is also necessary to rewrite the quadratic equation (j23H (whose roots are the sound waves 
A4 and A5) in terms of the wave frame quantities. Substitution of ([53]) and ([5"5"j) into ([231) yields 

< " ^ < = « 2 K , +K ?° tanhC • (69) 



Now we notice that according to Eq. (I10p we have 

Okq 1 d f w\ 1 1 d I ' w 



dp 2kq dp \ p J s 2 k' q cosh £ + n' n sinh £ d/5 



(70) 



where we have used Eq. (|27p by which the entropy is constant and thus w is only a function of 
p. We then substitute Eq. £70]) into (jSSJ) to find 



On - ° K + 

in which 



K., 



'2 



p d I w 

2dp\p 



«4 tanh£ = , (71) 



k' = yV 2 + w/p . (72) 

Equation (|7ip determines the phase velocity tanh £ in terms of the physical variables measured 
in the wave frame but since it is generally complicated depending on the explicit form of w{p) 
we can not go further. However, it is possible to continue for the ultra-relativistic case when 
KbT 3> mc 2 which implies that 

— =f-V , -=f-V , p=-w, 

w \T ) n a \T J 4 



by which it is easy to see 



w w _ 2 dP 1 W _ , , 

-rn" ' a= ^rw/' (73) 
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where the subscript "o" denotes the equilibrium point of the fluid at which it is at rest. By the 
above simplifications Eq. (fTTI) reduces to 

(«h - ^p){ k 'o + «4tanhf) = . 

Since n' n \ > 1 while | tanh£| < 1, the second factor in the above equation can not be zero 
and thus for the sound mode in the ultra-relativistic case we obtain 



w _ 



3 w _ 



(74) 



Here the upper (minus) sign indicates the case where the fluid velocity is negative with respect 
to the wave front which means that the wave runs faster than the fluid and thus it refers to the 
forward sound wave. Similarly the lower (plus) sign refers to the backward sound wave. 

As mentioned before there are only three independent equations namely Eq. (|68p when 
Eq. (|74p is substituted into it. It is more convenient to to write Eq. (|68p in the three orthogonal 
directions n, k'± and da./ dip. Thus, making the scalar product of (|68|) [after the substitution of 
((711) into it] by n yields 



i 



dp d£ 
T 



2~pl J I 2y 2% Jp( K 'l + %%.p)*P dip 



sinh£ 



■2 _i_ 3 wo z 
L + 2~%P 



. dn'± , dn 

--K ± • — h K ± • — , 

dip dip 



in which we have used the identity n'±_ • n = which gives 



dn' 



n + k 



dn 



dip dip 
Next, let us make the scalar product of ([68]) by k'±: 

, dn' i 



. 



dip y ^ ■■- L ' 2pl 

Finally the scalar product of (j68|) by dn/dip is 
dn' i dn 



si u I) i k';- + - T cosh £ A / \ >■"' 



dip dip 



dn 

dip 



sinhZ J k'I + -~p ^ cosh ^Jj^P ) l dn AM 2 



(75) 



(76) 



(77) 



(78) 



Thus we should solve the system of equations (75), ([771) and (fT8|) provided that n(i^) is a known 
function. 

If non of Eqs. (I77p and (|78p vanishes one can divide them by each other and after some simple 
vector calculations obtain 



dn' 



dn ( dn 



X — X K 



dip [dip \ dip 







which gives 



dn' 



dip 



Z(ip) x 



dn ( dn 



— x I — x K' 

dip \dip 



(79) 
(80) 
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where Z(</?) is an arbitrary continuous function. We will not go further in this way but alterna- 
tively seek fore more simple solutions. If we assume x = then it is possible to show after 
some calculations that this is not a consistent solution for the system of equations (75), (|77p 
and (|78|) . However, a consistent simple solution is found under the assumption 

dn' i , dn , , 

^r' n = - K -^ = °- < 81 > 

This condition with the help of (I77p gives 

Kj_ = const = R' ± . (82) 

By the use of ([81]) and ([82]) since | tanh£| < 1 we find a differential equation relating p to £ 
whose solution is 



-'2-2 

3w Q 



(cosh [±^=(£ - to) + cosh" 1 ( 1 + 

I L v3 V KjW 



1 , (83) 



where the upper (positive) sign refers to the forward and the lower (negative) sign denotes the 
backward sound wave. Now we should find k'±. It is clear from (|81|) and (|82|) that ^ x is 
perpendicular to both k'± and n thus ^ x is parallel to n x k'±. On the other hand the 
identity n • 4^ = and Eq. ([8"T|) imply that 4§ is also parallel to n x k'± . Therefore we conclude 



that 



dx,'± dn 

dip dip ' 



or equivalently 



k'±(¥>) = J\&)^^(kp' + K' ± (<p ) , (84) 

where 7r(y>') is an arbitrary nonzero continuous scalar function (we remember that ^ x can not 
vanish). Finally we substitute ([84"]) and ([82]) into ([7H]) to obtain 



^sinh^y k'I + =F cosh ^^Jj^P J = ~ 7T ( L P) • ( 85 ) 

Equations (I83p and (j85H are used to express both p and ^ as functions of 99 and this means that 
the problem is formally solved. Substituting of all physical quantities obtained above into the 
Lorentz transformation (|54|) provides all the things in the laboratory frame. 



5 Symmetry analysis for the vortex mode equation 

Before starting this section, let us mention that from here on we change all previous notations to 
quite new applications. So we forget the meaning of all letters or symbols used in all preceding 
sections and introduce new applications of them. 

We consider Eq. (|28p as a first order ODE and rewrite it in the following form 

dk ( k ■ n , \ , . 

n - 772-^T k )=0, (86) 



dt \ k 2 + w 

where w is a constant, t is treated as the wave phase, and k = (k%, fo, £3) and n = (n%,n2, n^) are 
some vectors in R 3 having the physical meaning of k and unit normal vector to the wave front 
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resp. We concern with the latter equation to find its point and contact symmetry properties 
and also give its fundamental invariants and a form of general solutions. 

It is notable here that in the mathematical structure of the simple wave solution, it is necessary 
to take the unit length for n (see Eq. (I13p ), However since Eq. (|86p is homogeneous and linear 
with respect to n, this condition is not essential in obtaining any solution. This condition 
appears important only for the compatibility of the simple wave structure. Regarding this fact, 
we make our symmetry analysis in both cases of arbitrary length and unit length for n and 
compare the results with each other. 

Throughout this section we assume that indices i,j varies between 1 and 3. Also each index 
of a function implies the derivation of the function with respect to it, unless specially stated 
otherwise. 



5.1 The point Symmetry of the Equation 



To find the symmetry group of Eq. (J86J) by Lie infinitesimal method, we follow the method 
presented in [32J. We find infinitesimal generators as well as the Lie algebra structure of the 
symmetry group of that equation. In this subsection, we are concerned with the action of the 
point transformation group. 

The equation is a relation along with the variables of 1— jet space J 1 (M,M 6 ) with (local) 
coordinate (t, k, n, q, p) = (t, ki,rij,q r ,p s ) (for 1 < i, j, r, s < 3), where this coordinate involving 
an independent variable t and 6 dependent variables fcj,n,- and their first derivatives q r ,p s with 
respect to t resp. 

Let M be the total space of independent and dependent variables. The solution space of 
Eq. (f86l) . (if it exists) is a subvariety 5a C J 1 (M,M 6 ) of the first order jet bundle of one- 
dimensional submanifolds of A4, that is, graph of functions ki,rij, of elements (t,ki(t),rij{t)) 

,, , . „ ™ , , ,. dki dk 2 dk 3 dm dn 2 

satisfying Eq. ([860 and the relations q\ = — , q 2 = <?3 = Pi = P2 = 

Ps = — — - are all fulfilled. 

Fi dt 

We define a point transformation on Ai with relations 

t = <f>(t,ki,rij), k r = Xr(t,ki,n,j), h s = tp s {t,ki,nj). 
where <f>, Xr and ip s are arbitrary smooth functions. Let 

i=l 

be the general form of infinitesimal generators that signify the Lie algebra g of the symmetry 
group G of Eq. (|86p . In this relation, T, Ki and Nj are smooth functions of variables t, k{ and 
rij. The first order prolongation [32j of v is as follows 

i j 

where K\ = DtQ\ + T q^t and Nj = DtQ\ + Tpjj, in which Dt is the total derivative and 

Q\ = Ki — T qi and Q\ = Nj — T pj are characteristics of vector field v [32]. By applying 
on (|86|). we obtain the following relation 

{ ((k 2 + to)(k • n + m) ~ 2 fcj) + Ni(k 2 + w)(k 2 +w- k 2 )] qi + [k h - (.Qj K ikj 

i j 

+pj K in .)- qi (r t - ^2 (qj T kj + Pj T nj y\ (n;(k 2 + w) - h(k ■ n)J | = 0, 
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whenever Eq. (|86[) is satisfied. We may prescribe t,ki,nj,q r ,p s (1 < i,j,r,s < 3) arbitrarily 
while functions T, Ki and Nj only depend on t, ki,rij. Thus Eq. (jB8"|) will be satisfied if and only 
if we have the following equations: 



^2 (»i(k 2 + w) - ki(k ■ n)j = 0, (89) 

i 

N^k 2 + w)(k 2 + w- kf) - ( n j( k2 + w ) - M k • n )) #j fc« 

^(k 2 + w) - A*(k- n)) T t + i^((k 2 + w)(k • n + m) - 2k^j = 0, (90) 

T ki (nj (k 2 + w) - kj(k • n) = 0, (91) 

T ki (nj(k 2 + w)- kj{k ■ nj + T k . (n^k 2 + w) - h{k ■ nj = 0, (92) 

( n j ( k2 + kj (k • n) Kj ni = 0. (93) 



4fe 



These equations are called the determining equations. From Eqs. (|90l) for each i we have 

Ni = (k 2 + w- kf)- 1 ^ (2 h t {k 2 + w)- (kn + m)\ Ki + ^ (n^k 2 + u») - %(k • n)) Kj 

j 

+ (n. i (k 2 +u;)-^(k-n))r t }. (94) 

Since without loss of generality, one may assume that ni 7^ 0. Also since k 2 + w 7^ so 

by Eqs. (f9Tj) and ([92]) we conclude that T just depends on t: 

T = T(t). (95) 

By solving Eq. (|89p with respect to i, we deduce the following relation of Ki s 

(^(k 2 + w)-k l (k-n)) K i = 0. (96) 



After differentiating of the latter equation with respect to rij when we apply Eqs. (|93h we lead 
to the following relations 

(k 2 + w - k 2 ) Ki-^ki kj Kj = 0. 

These relations suggest the general forms of K\ , K2 and K% as follows 

K 1 = K 2 = K 3 = 0. (97) 

By applying (j97H on relations (|94|) for different values of i, the forms of Ni s are also achieved: 

iV, = (k 2 + w-k 2 )- l ( ni (k 2 + w)-k l (k-n))T t . (98) 

Finally, the general form of infinitesimal generators as elements of point symmetry algebra of 
Eq. (I86p . which we call them point infinitesimal generators, for arbitrary functions T is as follows 

v = v T := T — + J2 {(k 2 + w - k 2 )-\ ni {k 2 + w) - h{k ■ n)) . (99) 

i=l 1 

The Lie bracket (commutator) of every two vector fields in the form of (|99|) straightforwardly is 
an infinitesimal operator in the same form of them. More explicitly, the commutator of operators 
vt and is vector field v T rp _ T y. Hence, the Lie algebra q = (vt) of point symmetry group 
G, when T is an arbitrary smooth function which depends on t, is a Lie algebra. 
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Theorem 1. The set of all point infinitesimal generators in the form of &9y\) is an infinite- 
dimensional Lie algebra of equation {86\) for arbitrary n (not necessarily unit). 

According to theorem 2.74 of [32J, the invariants u = I(t, ki,k2, &3, n\, n^, n^) of one-parameter 
group with infinitesimal generators in the form of (|99p satisfy the linear homogeneous partial 
differential equations of first order: 

v[I] = 0. 

The solutions of the latter, are found by the method of characteristics (See [32] and [33] for 
details). So we can replace the above equation by the following characteristic system of ordinary 
differential equations 

dt dk\ dk.2 dks dn\ dri2 dn 3 

By solving Eqs. (jlOOH of the differential generator (j99H . we (locally) find the following general 
solutions 

I a (t, k, n) = k a = d a , for a = 1, 2, 3 

i r l (101) 

^+3(*,k,n) = — jn^k 2 + w - kf,) -%(k-n)| = d p+3 , for = 1,2,3. 

where d a and dp are constants. The functions I\,l2, • • • , I& form a complete set of functionally 
independent invariants of one-parameter group generated by (|99p (see [32]). 

Similar to the theorem of section 4.3.3 of [33], the derived invariants (llOip as independent 
first integrals of the characteristic system of the infinitesimal generator (|99p provide the general 
solution 

S(t, k,n) := /i(/i(t,k,n),/ 2 (i,k,n), • • • ,I 6 (t 5 k,n)), 

with an arbitrary function fi, which satisfies in the equation v\p] = 0. This theorem can be 
extended for each finite set of independent first integrals (invariants) of characteristic system 
provided with an infinitesimal generator. 

In the following, we give some examples provided with different selections of coefficients of 
Eq. (|99p to show the method explicitly. We assumed that each appeared coefficient of vector 
fields is nonzero. 



Example 1. If we assume that T = 1 , then the infinitesimal operator (|99p reduces to the 
following vector field V\ = and the group transformations (or flows) for the parameter s are 
expressible as (t, k, n) — ► (t + s, k, n) which form the (local) symmetry group of v\. 
The derived invariants in this case will be as follows 

I a = k a , I/3 +3 = npik 2 + w - kj) - kpCk ■ n), for a, (3 = 1, 2, 3. 

Therefore the general solution corresponding to v\ when ji is an arbitrary function, will be 

S(t, k, n) = Jk, n(k 2 + w - kj) - k(k • n)) . 
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Example 2. Let T = t, then the infinitesimal generator is 

d ^ (9 
v 2 = t — + ^(k 2 + w- kfy^niik 2 + w)- ki(k ■ n)) — , 

j=i i 
Then, the flows of v 2 for various values of parameter s are 

(t, ki,n,j) -» (te s , k i} (k 2 + w - fc|) _1 {(n.,(k 2 + w) - kj(k • n)) e s + (k • n - rij kj) 
Also, we have the below invariants 

I a = k a , 1/3+3 = ^( ra ,s( k2 + w - kj) - kp(k- nfj, for a, f3 = 1, 2, 3, 
whenever defined and the general solution of Eq. (|86p as 

S(t,k,n) = /i(k, ^(n(k 2 + w - k 2 ) - k(k • n 
where [i is an arbitrary function. 

Example 3. In the case T = e t the infinitesimal generator (|99p changes to 

d ,„ 9 , 9n_1 , /„ 9 \ nn d 



v * = e * {qI + E( k ' + w - ^) _1 M k2 + «0 - • n)) — }, 



i=i " ' " J 

with group transformations of the parameter s transforming (t,ki,rij) to 

Jf-i „J\-i\ r„. i„ .„Jn „„t\-l(„ ru „ „ 1,-1/1,2 , „.. i„2\-i 



P(s) = (hi {e*(l - s e*)" 1 } , fc, , kj s e* (1 - s e*) -1 (n^ - (k • n - n^Xk 2 + u> - ]fc£ 

wherever defined. Independent invariants are 

I a = k a , 1/3+3 = e~* (raa(k 2 + w - k 2 ) - kp(k ■ n)), for a, /5 = 1, 2, 3, 

and hence the general solution of (I86p with respect to infinitesimal operator v 3 is an arbitrary 
function of these invariants. Indeed, if u = f(t,ki,nj) be a solution of Eq. (|86p then so is 
u = /(P(s)) for each s. 

When we assume that n to be of unit length, then by the action of «W (the first prolongation 
of the general form (|87p of infinitesimal generator v) on relation n\ + n 2 . + n 2 = 1 we tend to 
the following equation 

m JVi +n 2 N 2 + n 3 N 3 = 0. 



The last equation along with the deduced form of Ni s in (|98|) imply that Tt = 0. Hence T = c 
for arbitrary constant c and for each i, Ni = 0. Therefore, the form of infinitesimal generators 
reduces from relation (|99p to the below expression 

d 

V= dl- 

Theorem 2. The point Lie algebra of equation 186\) when n is a unit normal vector to the 
wave front is q = (Jj) isomorphic to the Lie algebra R. Therefore the point symmetry group of 
the equation with this additional condition is the group of phase translations. 



17 



5.2 The Contact Symmetry of the Equation 

In continuation, we change the group action and find symmetry group and invariants of Eq. (|86p 
up to the contact transformation groups. According to Backlund theorem [32J, if the number of 
dependent variables be greater than one (like our problem), then each contact transformation 
is the prolongation of a point transformation. In this subsection, we earn the structure of 
infinitesimal generators of contact transformations. The normal vector to wave front, n, is 
assumed to be either arbitrary (not necessarily unit) or of unit length and then we find the 
contact symmetry properties of the Eq. (I86p in these two situations. 

We suppose that the general form of a contact transformation be as following 

t = <p(t,ki,nj,q r ,p s ), ki = Xl{t,ki,n>j,Qr,Ps), n m = ip m (t,ki,rij,q r ,p s ), 

Qn — f]n (tj Tlj , q r , p s ) , p u — Q u (t,ki, Tlj , q r , p s ) , 

where i,j,l,m,n and u varies between 1 and 6; and (ftyXh^PrrnVn and £ u are arbitrary smooth 
functions. In this case of group action, an infinitesimal generator which is a vector field in 
J 1 (M, M 6 ), has the following general form 

B! -T! + E«£ + «A +ft £ + fl £}, (102) 

i=l 

for arbitrary smooth functions T, Ki, N m , Q m , P u (/ = 1, 2 and 1 < m, n, u < 3). 

Since our computations are done in 1-jet space, so we do not need to lift v to higher jet 
spaces and hence we act v (itself) on the Eq. (|86p . then we find the following relation 

3 

{Ni qi(k 2 +w)(k 2 +w- kj) + Qi(k 2 + w)( ni (k 2 + w) - fc;(k • n)) 

8=1 

-Ki qi (ni(k 2 + w) + (k • n)(k 2 + w - 2 fcj))} = 0. 

Since n / 0, so without less of generality, we can suppose that ri\ ^ 0, then the solution to this 
equation for would be 

3 

Qi = (k 2 +w)- 1 (n 1 (k 2 +w)-k 1 {k-n))- 1 {j2 ( K i ^{k 2 + w) + (k • n)(k 2 + w - 2 ki)) 

i=l 

3 

-Ni qi {k 2 + w )(k 2 + w- k 2 )) - Y, ft(k 2 + w)(ni(k 2 + w) - h(k ■ n))}. 

Therefore, the infinitesimal generator which we call it as contact infinitesimal generator is in the 
following form 

d ^ , d n l (k 2 + w) + (k-n)(k 2 + w-2k i )) 
V dt i (k 2 + w)(n 1 (k?+w)-k 1 (k-n)) Qi dq x 

EN ( ® k 2 + w - kf)) 9 \ p ^ (103) 

+ E Q 

3=2,3 



V dni m (k 2 + w) — k\ (k • n) dq± ) ^-r 1 dpi 
d nj{k 2 + w) + kj (k • n) d 



dqj n\ (k 2 + w) — k±(k ■ n) dqi J 
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Table 5.2 



The commutators table provided by contact symmetry. 





Vt 




Vis 


Vm 


Vn 


VP 


Vt 







VT + V Ka 


VT + V N/3 


VT + V Ql 


VT + V P 




— vt — 







VK a + V N0 


VK a + VQ 


VK a + V Prl 




—vt — 


VN 


-v Ka ~ v Nf3 





VN + V Ql 


v Nf3 + v Pll 


v Ql 


—vt — 


v Ql 


-VK a ~ V Q-y 


~VN fj ~ V Ql 





VQ-, + V Pv 




—vt — 




-v Ka - vp v 


-VNp ~ VQ V 


-v Ql - v Pv 






One may divide the latter form to the following vector fields, to consist a basis for Lie algebra 
q = (v) of contact symmetry group G: 

= d_ =K (— ^{k 2 + w) + {\L-n){\? + w-2ki)) _8_ 

VT &t' VK ' i \dk i + {k 2 +w)(n 1 (k 2 + w) -fci(k-n)) qi d qi 

„ d AT / d k 2 + w — k 2 )) d \ , 

V * =Pi Wi> ^ =Ar -fe- n l( k 2+W )- fcl (k.n) 9 ^)' (1 ° 4) 

_ / d rij(k 2 + w) + kj(k ■ n) d 
VQi = Qj \~dlfj ~ rn(k 2 + w) -ifei(k-n) Oil 

where 1 < % < 3 and 2 < j < 3. The commutator of every two of vector fields (|1U4|) is a linear 
combination of two operators (]104p which are generally in the form of those two operators again. 
Thus these vector fields construct a basis for the Lie algebra q of the contact symmetry group 
G. The commutator table is given in Table 15.21 for 1 < a, (3, r/ < 3 and 2 < 7 < 3. In this table, 
when the commutator of two vector fields is generally in the same form of some vector fields in 
(|104p . then we have used those general forms again to show the results of commutators. 



Theorem 3. The contact symmetry group of equation 186\) . is an infinite-dimensional Lie 
algebra generated by the contact infinitesimal operators jl04\ ) with the commutator table \5.2l 

One may repeat the above process for the problem of finding contact Lie algebra of Eq. (|86p 
with the supplementary condition of n to be unit. In this case there is another condition 
V [Eq. (|86p ] = 2(ni N\ + ri2 -/V2 + ri3 N$) = by the action of a contact infinitesimal generator on 
Eq. ([SB]) which must be added to other relations. Since n ^ 0, so we can suppose n\ 7^ and 
JVi = -2a N 2 - 2a iV 3 . Finally we can say 

Theorem 4. The contact symmetry group of equation (ESP consisting of unit normal wave 
front, is an infinite-dimensional Lie algebra and its Lie algebra is generated by the contact 
infinitesimal operators $103\) when we replace the coefficient N\ by —^N 2 — N3. The com- 
mutator table of these vector fields is in the form of Table \5.'£\ when we eliminate the row and 
column corresponding to UjVi and then change v^ 2 and v^ 3 to the new forms. 
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6 Summary and conclusions 



Towards a deeper understanding of the mysterious behavior of hydrodynamical equations it is 
necessary to look for various exact solutions. Among these solutions simple waves and multi- 
waves are the best for compressible flows up to present. These solutions show explicitly as a 
special case that how it is possible that smooth initial conditions convert to some discontinuities 
and singularities in future times. Thus, there is a significant hope that by a detailed and deep 
analysis of these waves one may find more general statements about the appearance of any 
non-smoothness regarding smooth initial conditions. 

In the present work a multidimensional version of simple waves introduced in References 
|15j and [T7] were employed for fully relativistic fluids and plasmas. Each wave front is a plane 
traveling with its own phase velocity vector. The intersection of different wave fronts is forbidden 
in the domain of the solution. Also at each instant of time there is a surface as the boundary 
between the two regions, the region of the validity of the solution and the forbidden region where 
the solution does not exist. This boundary generally moves and changes in the course of time. 

Similar to the nonrelativistic case[17J three essential modes were found, namely vortex, en- 
tropy and sound modes. Each mode suggests a wide variety of solutions while only very simple 
typical solutions were presented as some illustrations of the method of solving. Vortex and 
entropy modes were solved both in the laboratory and the wave frame. But, due to the high 
complexity of sound modes we studied them only in the wave frame. Further, as a special phys- 
ically valid example we considered the thermodynamically state equation at ultra-relativistic 
temperatures and obtained a complete formal solution in the wave frame. 

A symmetry analysis for the vortex mode equation (as a typical equation) led to the finding 
the structure of point and contact infinitesimal generators as well as fundamental invariants of 
the equation. In addition a form of general solutions implied by these invariants was obtained. 
Also we presented some examples for the point transformation case which tend to a precise 
determination of related symmetry groups. In the special case of our problem, the contact and 
point symmetry group of the vortex mode equation were both found to be infinite-dimensional 
Lie groups when the normal vector to wave front is not necessarily unit. When it is of unit 
length we find a one-dimensional point symmetry group while the contact symmetry group is 
still infinite-dimensional. The same procedure can be probably made for equations of other 
modes, namely the entropy mode and the sound mode. 
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